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^ ' Abstract. An integro-difFerential equation of hyperbolic type, with mixed 

, boundary conditions, is considered. A continuous space-time finite element 

method of degree one is formulated. A posteriori error representations based 
on space-time cells is presented such that it can be used for adaptive strategies 
based on dual weighted residual methods. A posteriori error estimates based 
on weighted global projections and local projections arc also proved. 
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1. Introduction 

We study the initial-boundary value problem (we use '•' to denote the time 
derivative) , 

u{x, t)-V ■ a{u; x, t) = /(x, t) in (0, T), 



>■ u{x,t)^0 onrDx(0,r), 

^ a{u; x^t) ■ n ~ g{x,t) onrNx(0,T), 

u{x,0) =u°{x), u{x,0)=v°{x) in rj, 



ly-^ ' that is a hyperbolic type integro-differential equation arising, e.g., in modeling 

, dynamic fractional/linear viscoelasticity. Here u is the displacement vector, / and 

CN ' g represent, respectively, the volume and surface loads. The stress a = a{u]x,t) is 

determined by 

a{t)=ao{t)- I K.{t- s)aois)ds, 

W ' with 







ao{t) = 2^oe(i) + AoTr(e(t))I, 

where I is the identity operator, e is the strain which is defined by e = ^(Vu + 
(Vu)-^), and Aq > are elastic constants of Lame type. The kernel K, is con- 
sidered to be either smooth (exponential), or no worse than weakly singular, so in 
both cases with the properties that 

(1.2) /C>0, /C'<0, |1/C|U,(K+) = < 1. 

For examples of problems of this type see, e.g., and |17j and their references. 
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We note that, e.g., completely monotone functions, i.e., functions b g ii(0, oo)n 
C^(0, oo), such that 

{-l)''D^b{t) > 0, t€ (0, oo), fc = 0, 1, 2, 

satisfy (|1.2p . when ||&||li(r+) < 1. Hence, Mittag-Leffler type kernels, that are 
weakly singular and arise in fractional order viscoelasticity, satisfy (|1.2p . see, e.g.. 

Existence, uniqueness and regularity of solution of a problem of the form (|l.ip has 
been studied in jl8| . A posteriori analysis of temporal finite element approximation 
of a parabolic type problem and discontinuous Galerkin finite element approxima- 
tion of a quasi-static (it sa 0) linear viscoelasticity problem has been studied, re- 
spectively, in 12 and |21| . For analysis and numerical solution of integro-differential 
equations and related problems, from the extensive literature, we mention [3], 
[I4] . [17], and their references. 

Here we consider a continuous space-time finite element approximation of de- 
gree one, cG(l)cG(l), for problem A similar method has been applied to the 
second order hyperbolic problems in [S], where the main objective is goal-oriented 
adaptive discretization of the problems. We introduce a posteriori error represen- 
tations, that can be a basis for goal-oriented adaptive strategies based on dual 
weighted residual (DWR) method, see [5] and [6] for details and computational 
aspects on the DWR method. For examples of application of the DWR approach 
for problems in solid mechanics see [10], [20], and the references therein. 

To evaluate the a posteriori error representations we need information about the 
continuous dual solution. Such information has to be obtained either through a 
priori analysis in form of bounds for the dual solution in certain Sobolev norms 
or through computation by solving the dual problem numerically. In this context 
we consider the former case and provide information through a priori analysis. To 
this end, we present two a posteriori error estimates based on global and local pro- 
jections. A weighted global a posteriori error estimate is obtained, using global 
L2-projections, for which the main framework is adapted from [8], and as an ex- 
ample of its specific application we refer to [9]. For the second a posteriori error 
estimate, that is based on local projections, we refer to (5j. 

The model problem (jl.ip is of hyperbolic type, and therefore our analysis are 
based on special treatment of discretization of the wave equation perturbed with a 
memory term. Here we only present the theory. 

The memory term causes the growing amount of data that has to be stored and 
used in each time step. The most commonly used algorithms for this integration 
are based on Lubich convolution quadrature |12j for fractional order operators, see 
also [19j for an improved version. For examples of the application of this approch 
to overcome the problem with the growing amount of data, that has to be stored 
and used in time stepping methods, see [1], [2], [3], [13], and [15]. In particular 
see, e.g., [1] where sparse quadrature together with appropriate a posteriori error 
representaions and estimates have been used in adaptive strategies based on the 
DWR approach. Wc plan to address this issue together with numerical adaptation 
methods such as DWR approach in future work. However, we should note that this 
is not an issue for exponentially decaying memory kernels, in linear viscoelasticity, 
that are represented as a Prony series. In this case recurrence relationships can be 
derived which means recurrence formula are used for history updating, see |10j and 
[21] for more details. 
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The outline of this paper is as fonows. In §2 we define a weak form of and 
the corresponding dual (adjoint) problem. In §3 we formulate a continuous Galerkin 
method of degree one. Then in §4 we present a posteriori error representations, from 
which, in §5 we prove a weighted global a posteriori error estimate. A localized a 
posteriori error estimate is also proved in §6. We remark on special cases regarding 
the memory term and regularity of the convolution kernel, when it is needed. 



2. Weak formulation and stability 

We let ft C M'', d = 2,3, be a bounded polygonal domain with boundary F = 
Fd U Fn where Fd and Fn are disjoint and meas(FD) ^ 0. We introduce the 
function spaces H = Lii^tY, Hr^ = L2{V^Y, and V = {v G H^{VlY : v\r^= 0}. 
Wc denote the norms in H and Hp^ by || • || and || • Htni respectively. 

2.1. Weak formulation. We define a bilinear form (with the usual summation 
convention) 

(2.1) a{v,w) = I {2iio€^j{v)e^.j{w) + \oeu{v)ejj{w)) dx, v,w 

which is coercive on V, and wc equip V with the inner product a(-, •) and norm 
\\v\\y ^ a{v,v). Wc define Au = — V • cro(i*), which is a selfadjoint, positive definite, 
unbounded linear operator, with 'D{A) = H'^ (il)'^ OV , and we use the norms \\v\\s ~ 
\\A'/^v\\. 

We use a "velocity-displacement" formulation of (jl.ip which is obtained by intro- 
ducing a new velocity variable. Henceforth we use the new variables ui = u, U2 = ii 
and u = (ui,U2) the pair of vector valued functions. Then the variational form is 
to find ui{t), U2{t) e V such that iti(O) = v? , ^2(0) = , and 

{ui{t),vi) - {u2{t),Vl) = 0, 

(2.2) {u2{t),V2) + a{ui{t),V2) - IC{t - s)a{ui{s),V2) ds 

Jo 

= {f{t),v2) + {g{t), V2)r^, yv,,v2 e V,te{0,T). 
Now we define the bilinear and linear forms B -.U x V ^ M. and L : V ^ M by 




(2.3) + (Wi,«i)- (U2,t.l)}dt+ (ui(0),Vi(0)) + (W2(0),W2(0)), 

L{v)= [{f,V2) + ig,v2)r^}dt+{u°,viiO)) + {v°,V2iO)), 

where 

(2.4) V , , ; V , , 7, 

V = [v = (wi, W2) : V e ^2(0, T; H) x ^2(0, T; V), Vi right continuous in t}. 

We note that if we change the order of the time integrals in the convolution term 
as well as changing the role of the variables s, we have the second variant of the 



4 



F. SAEDPANAH 



bilinear form that is, 

B{u,v) = / \{u2,V2) + a{ui,V2) - / IC{s ~ t)a(ui{t),V2{s)) ds 

(2.5) ^0 1^ Jt ^ 

+ (mi,wi) - (U2,vi)}rfi + (wi(0),vi(0)) + (m2(0),W2(0)). 

We use this variant of B in §4 for a posteriori error analysis. 
The weak form p.2p can be writen as: find u eU such that, 

(2.6) B{u,v)^L{v), VueV. 

Here the definition of the velocity U2 = ui is enforced in the L2 sense, and the initial 
data are placed in the bilinear form in a weak sense. A variant is used in where 
the velocity has been enforced in the sense, without placing the initial data in 
the bilinear form. We also note that the initial data arc retained by the choice of 
the function space V, that consists of right continuous functions with respect to 
time. 

Our error analysis is based on the duality arguments, and therefore we formulate 
the dual form of (|2.6|) . To this end, we define the bilinear and linear forms B* : 

V* xU* ^ R, L*:V*^ R, for r e R^", by 

B*{v,z)=J {vi,zi) + a{vi,Z2) ~ IC{s - t)a{vi,Z2{s)) ds 

(2.7) - {V2,Z2) - {v2,Zi)}dt+{vi{T),Zi{T)) + {V2{T),Z2{T)), 
L;{v)= £{ {V, + {V2 ,j2)}dt+ {V, (T) , zf) + {V2 (T) ,z^), 

where ji, j2 and zf , represent, respectively, the load terms and the initial data 
of the dual (adjoint) problem. In case of t = 0, we use the notation B*,L* for 
short. Here 

(2.8) V . . y V , , 

' V* = {v = ivi,V2) ■■ V e i2(0,r;y) x L2iO,T;H),v, left continuous in t}. 

We note that, recaUing (HH), U C V*, U* C V. 

We also note that B* is the adjoint form of B. Indeed, integrating by parts with 
respect to time in B, then changing the order of integrals in the convolution term 
as well as changing the role of the variables s, t, we have, 

(2.9) B{u,v) = B*{u,v), yueU,veU*. 

Hence, the variational form of the dual problem is to find z € U* such that, 

(2.10) B*{v,z) ^ L*{v), VweV*, 
that is a weak formulation of 



T 



d 



(2.11) Z2+Az2~J JC{s~t)Az2{s)ds = ji- —j2, 

with initial data z'[,Z2, and function j = (ji,j2) that is defined by L*{w) 
Jo (w, j) dt. 

We note that, since is a solution of (|2.10l) . we have 

(2.12) Zi^-Z2-j2- 
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2.2. Stability of the solution of the dual problem. Wc know that stability 
estimates and the corresponding anlysis for dual problem is similar to the primal 
problem, however with opposite time direction. Hence, having a smooth or weakly 
singular kernel with (|1.2[) . we can obtain the stability estimates (|2.13p from, e.g., 
[4] or [11] for the continuous dual solution. Here we state stability estimates for the 
continuous dual solution in the following lemma, and we omit the proof for short. 
We note that the stability constant in ()2.13|) does not depend on t, and Gronwall's 
lemma has not been used, see H] and [TT]. See also (TB], [T7] and HI], where stability 
estimates have been represented, in which the stability factor depends on t, due to 
Gronwall's lemma. 

Lemma 1. Let z be the solution of the dual problem (|2.10[) with sufficiently smooth 
data zf , Z2 , ji, j2- Then, for some constant C ~ C(k), we have stability estimates 

(2.13) ||zi(t)||, + < C{\\zf\\i + + ^"^(lljill, + ||j2||m) dr}. 

We note that, if we set ji = .72 = in (|2?T3|l and recaU ([2TT|) - ((2?T2l) . then for 
r G [t, T] and I €R, we have 

= ||i;2W|l;-i = II - Az2{r) + /C(s - r)Az2{s) ds\\^_^ 

< II Az2(r)!|z-i + ^ IC{s - r)\\Az2is)\\i.i ds, 

that, recalling ||/C||li(r) = ^or some C — C{k) implies 
\\ziir)\\i-i ^ ||i;2(r)||;_i < (1 + k) max Pz2(r)||,_i < (1 + k) max ||z2(r)||/+i. 

From this and (I2.13P wc conclude the stability inequality 

(2.14) ||ii(t)||/-i + \\zi{t)\\i + ||z2(t)||/+i < C{\\zJ\\i + \\z^\\i+i}. 

3. The continuous Galerkin method 

Here we formulate the continuous Galerkin method of order one, cG(l)cG(l), 
for the primary and dual problems ()2.6p and ()2.10p . 

3.1. The cG method. Let = to < ti < • • • < tn-i < t„ < • • • < tw = T be a 
partition of the time interval [0,T]. To each discrete time level i„ we associate a 
triangulation 7^" of the polygonal domain with the mesh function, 

(3.1) hn{x) = hx = diam(A'), x G A', K e 7^", 

and a finite element space V^' consisting of continuous piecewise linear polynomials. 
For each time subintcrval /„ (tn-ij^n) of length = — t^^i, we define 
intermediate triangulaion 7^" which is composed of mutually finest meshes of the 
neighboring meshes 7^", 7^"~^ defined at discrete time levels t„, t„_i, respectively. 
The mesh function /i„ is then defined by 

(3.2) hnix) =hK = diam(7s:), x e K, K e T^- 

Gorrespondingly, we define the finite element spaces Vj^ consisting of continuous 
piecewise linear polynomials. This construction is used in order to allow continuity 
in time of the trial functions when the meshes change with time. Hence wc obtain a 
decomposition of each time slab fi" = ilx/„ into space- time cells K"" — Kxin, K £ 
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7^" (prisms, for example, in case of C M^). We note the difference between the 
mesh functions /i„ and and this is important in our a posteriori error analysis. 
The trial and test function spaces for the discrete form are, respectively: 

U continuous in il x [0,T],U{x,t)\j^ linear in t, 

V{-,t) continuous in n,V{-,t)\i^ € (V^)^, 
V{x,t)\i^ piecewise constant in i|. 

We note that global continuity of the trial functions in Uhk requires the use of 
^hanging nodes^ if the spatial mesh changes across a time level t„. In these 'irregular' 
nodal points the unknowns are eliminated by interpolating values at neighboring 
'regular' nodal points, see [5], [35] and the references therein for practical aspects. 
We allow one hanging node per edge or face. On the other hand, the test functions, 
which are allowed to be discontinuous in time, arc defined in the time slabs 11 x /„ 
on the spatial meshes corresponding to t„. In this discretization trial functions have 
support in both time intervals /„ and In+i adjacent to time instants Therefore 
on each time interval /„, the trial functions defined at time instant on a mesh 
TJj"^ overlap with test functions defined at time instant t„ on a mesh 7^" . Hence, 
to form space-time integrals of multiplication of trial functions and test functions, 
we need the 'union' of the two triangulations, that we have denoted by 7^". We also 
note that, due to the structure of the elements of 7^", computations are feasible with 
reasonable effort if the grids 7^"~^ and 7^" are related. Hierarchically structured 
spatial meshes, where 7^" is the set of most refined cells from the two grids, reduces 
the substantial work caused by the mesh transfer from one time level to the next. 
More details and references on the practical implementation can be found in 'S' ■ 

In the construction of lAhk and Vhk we have associated the triangulation 7^" with 
discrete time levels instead of the time slabs fi", and in the interior of time slabs we 
let U be from the union of the finite element spaces defined on the triangulations at 
the two adjacent time levels. Associating triangulation with time slabs instead of 
time levels would yield a variant scheme which includes jump terms due to discon- 
tinuity at discrete time leveles, when coarsening happens. This means that there 
are extra degrees of freedom that one might use suitable projections for transfering 
solution at the time levels see [TT| . 

The continuous Galcrkin method, based on the variational formulation (|2.2p . is 
to find U G Uhk such that, 

(3.4) B{U,V)^L{V), -iVeVuk- 

The Galerkin orthogonality, with u = (ui,U2) being the exact solution of ()2.2|) . is 
then. 



Uhk = [u^{Ui,U2) 



(3.3) 



Vhk = [v = {VuV2) 



(3.5) 



B{U-u,V)^Q, yv^Vhk- 
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From p.4|) we can recover the time stepping scheme, 
{{Ui,Vi) - {U2,Vi)} dt ^ 0, 

^ {{U2,V2) + a{Ui,V2)- J^IC{t-s)a{Uii3),V2{t))ds}dt 
{{f,V2)dt+{g,V2)r^}dt, yiVi,V2) eVhk, 



(3.6) 



for suitable choice of u° , u° € as approximations of the initial data w*^, v". 

We define the orthogonal projections Vh^i : H ^ V^^ and Vk,n ■ L2{InY ^ 
Po(/„), respectively, by 

(3.7) 



with Pq denoting the set of all vector-valued constant polynomials. Correspond- 
ingly, we define VhV and VkV for tGl„{n = 1, - ■ ■ , N), by {Vhv){t) = Vh,nv{t) and 
T'kV = Pk,n{v\i„). We note that, as a natural choice, we can set the initial data in 
(13^ as 



(3.8) ul = Vhu°, = VhiP. 

We also note that, when we do not change the spatial mesh or just refine the 
spatial mesh from one time level to the next one, i.e., V^~^ C V^, n = 1, ... ,7V, 
then we have = V^- 

We introduce the linear operator An^r '-^h ^"7 

We set An = A„^„, with discrete norms 

\\vn\\h,i = ll^J/^Wnll = \/(w„,A>„), «n € V;^ and Z € R , 
and Ah so that A/iW = AnV for w € V^. We use Ah when it acts on . 

4. A posteriori error estimation: error representation 

Having certain regularity on the data, i.e., initial data vP ^ and the force terms 
/, (7, there are still two types of limitation for higher global regularity of a weak solu- 
tion of . One is due to the mixed Dirichlet- Neumann boundary condition. This 
type of boundary condition are natural in practice, and a pure Dirichlet boundary 
condition cannot be realistic in applications. Other limitation is the singularity of 
the convolution kernel /C. This means that even with the pure Dirichlet boundary 
condition, higher regularity of a weak solution is limited, see [11], [18], though with 
smoother kernels we can get higher regularity. These, and other general motiva- 
tions such as no practical use of a priori error estimates, call for adaptive meshes 
based on a posteriori error analysis. 

Here a space-time cellwise error representation is given. The main framework is 
adapted from [6], and a general linear goal functional i*(-) is used. For an example 
of a global goal functional, see ExamplelH This error representation can be used for 
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goal-oriented adaptive strategies based on dual weighted residual (DWR) method. 
For more details on dual weighted residual method and its practical aspects for 
differential equations, see [S], [B] and references therein. 

In the following, we recall and define some notations for this section. Wc denote 
the space-time cells i^" = A' x /„ and dK"^ ~ dK x /„, for X e 7^", and it 
should be noticed that dK^ is not the boundary of if". Also for a simplex A', we 
define the inner products {■,-)k^ = Jj {■,-)k dt, {■,-)dK" — Jj {■,-)dK dt, and the 
corresponding A2-norms 

Mk" = (y^ \\-\\ldt)'^ =[j^ \\-\\l,^K)dt)'\ 

II • hK^ = (/ II • 11^^ '^^T' = (/ II • ll'^(e^) '^^T'- 

We also denote 

(4.1) {IC*vy{t)= IC{t - s)v{s) ds. 

Throughout we use the usual notations a A 6 = min{a, b} and ay h = max{a, 6}. 

Now we present three a posteriori error representations, and we note that the 
second error representation. (|4.3p . is space-time cellwise. 



Theorem 1. Let u and U be the solutions of (j2.6|) and (j3.4p . respectively, and 
L*(-) be the linear functional defined in (|2.7p . Then, denoting the error e ^ U — u, 
we have the error representations 

N 4 N n 

(4.2) L*ie)= y: qo.k+y: e E0"^+ee e 

N 5 

(4.3) L*{e)= 0o,A- + E E J2®Ik, 

AT 4 N N 

(4.4) A*(e)= Y eo,K + E E E0-- + EE E 

where, with z^k € Vhk being an approximation of the dual solution z and EhkZ = 
Zhk — z being the error operator, 

eo,K = (c/i(o) -u",£;,,fezi(0))^, + (c/2(o) -f°,s,,fcZ2(o))^, 

= (Ui - C/2, S^fcZi)^f , = (f/2 - /, EhkZ2)K- , 

©3,^ ~ 9, EhkZ2)dK^, 04,K = {'Td, EhkZ2)dK", 

(4.5) ©s:^ = ((^ * rrf)^ (0, A,,Z2(t))^^ dt, 
e5K = -(rdJ JC{s~t)EhkZ2{s)ds) 

ef^: = - / irS), / ' '»C(s - t)S„fcZ2(s) ds) dt. 
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Here rd are the residuals representing the jumps of the normal derivatives ao{Ui)-n, 
which are determined by 



(4.6) rrfir 



-i[(7o(t/i)-n] tfTcdK\dn, 

«/ r c dn, 



and Qd are the contribution from the Neumann boundary condition defined by 
(4.7) 

r([/i) • n = (ao(C/i) - /o/C(t - s)ao{Ui{s)) ds) ■ n if V d dK n T^,, 



9d\T 



^ otherwise. 

Proof. Using the identity (j2.9|) and the Galerkin orthogonality p.5|) we have, 

L*ie) = B*ie,z) = B{e,z)^B{e,Ehkz) = B{U,Ehkz)-B{u,Ehkz) 
= B{U, Ehkz) - L{Ehkz) = R{U; E^kz), 

where R(U; •) is the residual of the Galerkin approximation J7 as a functional on 
the solution space U*. Then, by the definition of B and L in (|2.3p . we have 



(4.9) 



L*{e) = (C/i(0)-M°,S,,fc2i(0)) + {U2{0)-v'',EhkZ2{0)) 

+ / {{ill , EhkZi) — {U2, EhkZi) + ({72, EhkZ2) + a{Ui, EhkZ2) 

/C(t - s)a(C/i(s), £;,,feZ2(t)) dsj 
|(/, £'/ifeZ2) + (.g, -Bhfc22)rN} c^*- 



Now, by partial integration with respect to the space variable and recalling rd from 
(14.61). we obtain 



(4.10) 

a{Ui,EhkZ2)dt 



N „ N 



/ a{Ui,EhkZ2)Kdt^^^ I ^ [ao{Ui) ■ n,E,ikZ2)gj^dt 

n=l "'-f" _R-er^.. n=l "'■f" Kef;" 

N „ 

{J2 {-[MUi)-n],EhkZ2)E+ J2 {MUi) ■ n, EhkZ2) e} dt 

N „ 

X! / X! {i'^d,EhkZ2)dK + {(ToiUi)-n,EhkZ2)gj^^j.^^dt, 



where f^, respectively, the sets of the interior edges and the edges on the 

Neumann boundary, corresponding to the triangulation 7^. 
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For the convolution term in ()4.9p . similar to (|4.10p . we have 



i-T rt 

IC{t^ s)a{Ui{s),EhkZ2it)) ds dt 



'0 Jo 



IC{t-s)Y, a{Ui{s),EhkZ2{t))^dsdt 
E/ E/ ^(i-^) E {-Wo{Ui)-n],E,,Z2{t)) 



(cro(C/i) • n, £;,ifeZ2(t)j \ ds 



N „ n tjAt 



^ / /C(t-,s) 5] {(r,(.),i?,,Z2W 
+ (o-o(?7i(s)) • n, Z2(i)) [ ds di, 



that implies 



T /.t 

K.{t- s)a{Ui{s),EhkZ2{t)) ds dt 



"'0 



N . n t,At 

E/ {[ IC{t-s)rdis)ds,EnkZ2it) 



tjAt 



dK 



dt 



+ ( IC{t-s)ao{Uiis))-nds,EhkZ2it)) } 

EE E / (/ IC{t-s)ra{s) ds,EhkZ2it)) dt 



y2 [ ( [ ICit~s)cToiUi{s))-nds,EhkZ2it)) dt 



EE E / (/ /C(t-s)rd(s) ds,S^feZ2(i)) di 
E/ E {J^ ^t^s)c7oiUiis))-nds,E,,kZ2{t)^^^^^ dt. 



n=l J" ifer^ 
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This, together with (|4.10p and recahing gj, from (|4.7|) . imply 

N 



I] / (/ IC{t-s)rd{s) ds,EhkZ2] dt. 



Hence, from ()4.1ip and space-time cehwise representation of the other terms in 
(|4.9p . wc conclude the first a posteriori error representation (|4.2p . 

Now wc prove the second error representation (|4.3p . To this end, in (|4.8p we use 
p.Sp the second variant of the bilinear form B, in which the convolution term has 
been rearanged. Therefore we just need to study the rearanged convolution term 
that, similar to ()4.10p . is writcn as 

JC{s - t)a{Uiit),EhkZ2{s)) ds dt 

N 



(4.12) 



V / V a(c/i(t), / IC{s-t)EhkZ2{s) ds) dt 

/ y2 [rd, / IC{s - t)EhkZ2{s) ds) dt 
+ 11/ {MUi)-n, JCis-t)EhkZ2{s) ds) dt. 

n=l-^^" '"^ 



For the second term of the right side, we exchange the role of the variables s, t, and 
we change the order of the time integrals to obtain 

y2 (o-o(J7i) • n, / K,{s - t)EhkZ2{s) ds\ dt 

„=lJ/„^ -It ^Tn 

(4.13) = / / nt~s){aQ{Ui{s))-n,EhkZ2{t))^ dt ds 

Jo Js " 

- V / ( / JC{t - s)ao{Ui{s)) ■ n ds,EhkZ2{t)] dt, 
that, with (|4.12[) . we have 

T i-T 

lC{s - t)a{Ui{t), EhkZ2{s)) ds dt 



Jt 



+ f / /C(i-s)(To(i7i(s)) -n ds,£;ftfcZ2(0) | dt. 
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This, together with (|4.10p . imply 

[a{Ui{s),EhkZ2[t))~ j lC{s-t)a{Ui{t),EhkZ2{s)) ds} dt 

N . 

(4.14) = X! X! Wd,EhkZ2)dK^ + {9d,EhkZ2)dK^ 

- {^'d, IC{s - t)EhkZ2{s) dsj 

that is the replacement for (|4.1ip . 

Hence, from (i4l4|) and space-time cellwise representation of the other terms in 
()4.9p . we conclude the second error representation ()4.3p . 

Finally, we prove the third error representation (|4.4p . In (|4.8p we use (|2.5|) . the 
second variant of the bilinear form B, and similar to (|4.10p we have 

T 

/C(s — t)a{Ui{t), EhkZ2[s)) ds dt 

N „ N „ti 



= f 111' IC{s-t)a{Ui{t),EhkZ2{s)) ds dt 

n-l i-n-^'j-iVt 
4.15 ^ 

N „ N „t 

„=i -^-fr. j=„ -^tj-ivt A'ef^^ 

+ fTo(C/i(t)) •n,£;,,feZ2(s)) I 
V / SA'nrN J 

Then, rewriting the second term in the brace similar to ()4.13p . we have 

T i-T 

K.{s - t)a{Ui{t), EhkZ2{s)) ds dt 



dK 



ds dt. 



Jt 



N N 



EE E / / ICis-t){r,it),E^kZ2is)) gji ds dt 
E/ E (/ ICit-s)aoiUi{s))-nds,EhkZ2it) 



N 

" ^ ^ dt 

This, together with (|4.10p . imply 

' [a{Ui[s),EhkZ2{t))^ ICis-t)a{Ui{t),EhkZ2{s)) ds} dt 

N 

(4 16) "EE {(^d'-^'«fe^2)aA" + (5d,-E^/ife22)aA"} 

n=life7T> 

Af Af 



EE E / (r<i(t),/C(s-i)£;,,fc^2(s)) rfsdi. 



that is the replacement for (|4.1ip 
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Hence, from ()4.16|) and space-time cellwise representation of the other terms in 
()4.9p . we conclude the third error representation ()4.4|) . Now the proof is complete. 



We note that aU error indicators in (|4.5[) depend on the unknown dual solution 
z, which is not available. For strategies to evaluate the error indicators, that are 
based on suitable approximation of the dual solution, we refer to [5] and [6] . When 
the spatial meshes change from a time level to the next one, evaluation of functions 
is carried out by means of interpolation/extrapolation. We recall that one hanging 
node is allowed per edge or face. 

We also note that the difference between a posteriori error representations (|4.2p - 
(|4.4p are the error indicators 65 65 and 6^j^, where we apply the convolution 
integral either on the residual or on the error term EhkZ2- Obviously we choose 
the cheaper one, depending on the method we choose for computing or estimating 
EhkZ2- It should be noted that the first and the third error representations (|4.2p 
and (|4.4p are not space-time cellwise. 

5. A posteriori error estimates based on global projections 

In order to evaluate the a posteriori error representations (|4.2P " ()4.4p . we need 
information about the continuous dual solution z. Such information has to be ob- 
tained either through a priori analysis in form of bounds for z in certain Sobolev 
norms or through computation by solving the dual problem numerically. In this 
context we provide information through a priori analysis and we leave the investi- 
gation on the second case to a later work. 

Sometimes the target functional one is interested in is (almost) global. Then, 
instead of numerically approximating the dual solution, one can get cheaper error 
indicators by using analytical a priori estimates for the dual solution. We first 
present a weighted global a posteriori error estimate, using global L2-projections 
'Pkj'Ph defined in p.7p . and error estimates of Vh in a weighted L2-iiorm. 

We recall the weighted global error estimates of the L2-pi'ojection Vh p.7p . see 
[7]. First we recall some notation. Let 7~ be a given triangulation with mesh 
function /i, and for any simplex K £ T , pk denote the radius of the largest ball 
contained in the closure of K , that is K. A family F of triangulations T is called 
non-degenerate, if there exsists a constant cq such that 



Let Sk = {K' € T K' n K 0}. Also Sj- and (5jf be measures for a given 
triangulation T and a given family respectively, defined by 



□ 



cn = max max . 

TeJ^ KeT PK 



St ~ max max 11 — hi^, /hi, 
KeTK'eSK K / « 




We define the error operators Ehk, Eh, and Ek by 



EhkV = {VkVh - I)v, Ehv 



[Vh - I)v, Ekv 



{-Pk - I)v, 



and we note that 



(5.1) 



Ehk — Eh + EkVh- 
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We recall standard error estimation for the time projection error Ek, that is, 



(5.2) J \Ekv\ dt^ J \Ek,nv\ dt < Ckl J \d^v\ dt, 7 = 0, 1. 



We also quote the error estimates for the spatial projection error E/^ from [7], which 
is stated in the following lemma. 

Lemma 2. Assume that the family T of traingulations T he non- degenerate. Then 
for sufficiently small Sjr, there exists a constant C such that for any triangulation 
T € we have, 



(5.3) Wh-'EhvW <C\\V''v\\, s = 0,1,2, Vw e ff", 

(5.4) \\h-'VEhv\\<C\\V'+'vl s = 0,l, V^;Gi^^ 

where V ' denotes the usual gradient. 

For more details on the practical aspects of (5jr, see f^. 

For the next theorem we recall the mesh functions /i„ , hn from p.ip , p.2p , and 
we define the notations 



min Hk, h^ax.n = max Hk, 



(/ IC{s ~ t) dsy^^ , /C„j=(/"' IC{s~t)dsy^^ 



We also use the scaled trace inequality, for any simplex K G Th- 



(5.5) \\v\\0K < CihJ,'/^\\v\\K + h'f\\Vv\\K). 



Theorem 2. Let u he the solution of (|2.6p . and U he the solution of (|3.4p with 
a non- degenerate family Th of triangulations 7^", 71 = 0, 1, ... , N , with sufficiently 
small Sjr^, such that the weighted glohal error estimates (|5.3p and ()5.4p hold. We 
also assume that K, e Li(R"'"). Then, denoting the error e = U ~ u, we have the 
weighted a posteriori error estimate, for a = 0, 1, 2, /? = 1,2, 7 = 0, 1, 

|L*(e)| < C max^ {|lV"zi(t)|l, |lV'^Z2(t)|Miar'^i(i)IM197^2(t)|l} 
(5.6) N 



n=l ■'^^ 
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whe 



T, = \\K{Um-u')\\ + \\h',{Um-v')l 
Tn^\\K{Ui-U2)\\ + \\fii{U2- f)\l 

1/2 



Tfc = k'tWEk{iJ^ - c/2)|| + fc,;j(||SfcAc/iii 

KeT" 



with Cn{fi) = /lf„™,„, Cn,w(/3) = Kn£n^n,TY.f=n><^n.3'^'Lalj- 

Proof. Let z e V* be the solution of the dual problem p.lOp . From the definition 
of the L2 projeetions VkjVh in (\3.7\i and the test space Whk in p.3p we have 
VkVhZ e W/ife. Therefore, using and ([ET|) we have, 

(5.7) L*{e) = R{U- E^kz) = R{U- E^z) + R{U; EkVhz). 

We study the two terms at the right side of this equation. 

For the first term, recalling (|4.8p and using the a posteriori error representation 
(|4.3p . we have 

R{U;Ehz)^ L*{e) 

= E {(C^i(0)-M°,i?^zi(0))^ + (C/2(0)-«°,£;„z2(0))^} 

(5.8) +E/ E |(^i-^2,i^/.^i),, + (t/2-/,i?^.^2)i^ 

+ (r^, EhZ2)dK + (ffd - .g, EhZ2)dK 



/I N .£ U 

/C(s-i)£^,,02(s)rfs) Ut = ^I,+^II, 
■^^J 1=1 i=l 



We then, for each term, use the Cauchy-Schwarz inequality twice. First on the 
local elements K and dK, to obtain local L2-norms, and then on the sum over 
the elements to obtain global norms such that the weighted global error estimates 
(|5.3p - (|5.4p can be used. For Ii, using (|5.3p . we have, for a = 0, 1, 2, 

(5.9) Ii < ||/i^(C/i(0) - 1.0)1111 V"i?h^i(0)|| <C||;if^(C/i(0)-7/0)|||lV"zi (0)11, 
and in a similar way, for /3 = 1,2, 

(5.10) l2<qi/io^(C/2(0)-^;°)|l||V^Z2(0)|l. 
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Now we estimate the second terms 11;, i = 1,...,5. Using the Cauchy-Schwarz 
inequaUty and the error estimate (|5.3|) . we have, for a = 0, 1, 2, 



III = ^ / ~ U2,EhZ,)^dt 

(5.11) <^ / \K{u,-U2mKrEHZi\\dt 

1 'J In 



N 

<Crnax||V"zi(t)||^ / \\KiUi - U2)\\ dt, 

' n— 1 



and similarly we obtain, for /3 = 1, 2, 

N 

(5.12) II^^CCmaxllV^z^WllE / \\h?.{U2 - DW dt 

For II3, we first have, 

N 



Then, by the scaled trace inequality (|5.5p and the weighted global error estimates 
(I531)-(l5lll), we obtain, for /3 = 1, 2, 

E h-K\\EhZ2\\lK < C E {^7/P'^^2||?, + /i;,'||Vi?,Z2||?,} 

These imply the estimate, for /? = 1, 2, 

(5.13) Il3<Cmax||V^Z2(0llEy/™'n( E ^i'lkdllk) dt. 



In a similar way, we have, for /3 = 1,2, 



(5.14) Il4<Cmax|lV^Z2(t)|jE / hit,n{ E ^^Il5<i - sllk") c?*- 
Finally we study II5. To this end, first we note that, 

IdK 

Ken 



"^^E / ( E /^i'lk<i(t)| 

E J IC{s~t)EhZ2{s)ds 



2 s 1/2 

dt. 

dK 
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Then, using the Cuachy-Schwarz inequahty, we have 



j-T 

h-^^ J IC{s-t)EhZ2{s)ds 



OK 



< Y ~^k' [ IC{s-t)ds [ lC{s - t)\\EhZ2{s)\\lj, ds 
Ken- 

<K.It j^lC{s~t) Y hK'\\EHZ2i.s)\\li,ds, 



that using the scaled trace inequahty (|5.5p and (|5.3p - (|5.4p . we have, for /? = 1, 2, 



2 

dK 



Y ^^i/ 'C{s-t)EhZ2{s)ds 

\\V EhZ2{s)\\\] ds 

+ \\h-if+^\/EhZ2{s)\\\]ds 

N .f. 

< CKLn'^lrY'^Z:!^ / ^(.^ - t)l|V^Z2(s)f ds 



IC{s — t) ds. 



Hence we have 



II5 < Cmax||V^Z2(0ll 

[0,T] 



(5.15) 



N 



71=1 



N 



E L {h;;^'ni^n.Tj2K^,,i^n.,){ E h'M\lKy^'}dt. 
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Putting ([5J|) - (|5?T5l) in ([SH) we conclude, for a = 0, 1, 2, ^ = 1, 2, 

R{U;Ehz)<Cm^{\\S/"zi{t)l\\V^Z2{t)\\} 



(5.16) 



x\\\h^{um-^n\\ + \\h^o{U2iO)-v'')\\ 

N „ 

E / {fKiUi-u,)\\ + mu2-f)\\ 

n=l 



1 1/2 ^ 

Now we study the second term R{U]EkPhz) in (|5.7p . We have 

R(f/;^feP,.z) = R(;7; ± / {{Vkf . EkVhZ2) + [Vug, EkVhZ2)r^] dt 

Jo 

= {Ui{0)-u'',EkV,MO)) + {U2{0)-v",EkVhZ2{0)) 

N „ 

+ E / {{Ui-U2,EkVhZi) +a{Ui,EkVhZ2) 



n—1 



JC{t - s)a{Ui{s), EkVhZ2) ds + {U2 - Vkf, EkVhZ2) 
- {rkg,EkrhZ2)T^ + {Ekf,EkVhZ2) + {Ekg,EkVhZ2)r^]dt. 

Recalling the initial data p.8l) . C/i(0) = 'PhUi{0), i = 1,2, the first two terms on the 
right side vanish. Besides, from the second equation of p.6p we have, for V S Vj^, 

{ {U2, V) - {Vkf, V) - [Vkg, V)r^ } dt 



[a{VkUi,V)^a{Vk K.{t - s)Ui{s) ds,v)] 



dt. 



Hence, we conclude 



(5.17) 



R(C7; EkVhz) = [ { (^'^^^i " ^2), EkVhZi) - a{EkUi, E„VhZ2) 

n=l 

+ a[Ek j^]C{t ^ s)Ui{s)ds,EkVhZ2 
+ {Ekf, EuVhZ2) + {Ekg, EkVhZ2)v^ } dt. 
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For the last term, we have 

TV 



J2 I iEkg,EkVhZ2)r^dt 

n=l ■'I^ 

^ E / ( E ^A-'ll^fe.9llk-) ( E hKWEkVH 

n=l -f" R-cT" R'pT" 



By the scaled trace inequality ()5.5p and local inverse inequality, 
we have 

J2 hK\\EkVHZ2\\lKnr.<C ^ {|li?.n^2||i- + /4llVi?.n^2|||} 

<C ^ {||i?.P.Z2|li- + llM.Z2|li^}=C||£;feP,z2f. 

Hence 

W „ ^ /■ 1/2 

^ / (i?fe5,i^fc7'„z2Kdt<C^ / ( J2 ^KWEkgWlK) \\EkVhZ2\\. 
Considering this in (j5.17p . using the Cauchy-Schwarz inequality we have 

RiU;EkVhz)<cY, [ \\\EkiUi-U2)\\\\EkVhZi\\ + \\EkAhUi\\\\EkVhZ2\\ 

+ \\Ek f IC{t~s)AhUi{s)ds \\EkrhZ2\\ + \\Ekf\\\\EkVhZ2\\ 
II Jo 

+ ( E h],'\\Ekg\\lKy^'\\E,rhZ2\\\. 
iter," ^ 

From this, together with L2-stability of the L2-pi'ojection Vh and the error estimate 
((5?2)) for Eu^Vk- I, we conclude, for a = 0, 1, 2, 7 = 0, 1, 

N 

R{u-EkVnz)<c^^^{\\dr'z,mA\diz^m]Y. / {fcr'iif>i-f^2ii 



(5.18) +kl\\EkAhUi\\+kl\\Ek / /C(t-s)i,X/i(s)ds| 

1/2^ 



+ kl\\Ekf\\+kl{^ ^ } 



if ST" 



Hence, putting (|5.16p and (|5.18p in (|5.7p . we conclude the a posteriori error 
estimate (|5.6p . and this completes the proof. □ 



Wc note that, to obtaining a computable error bound in (|5.6p . one aspect is 
to estimate or eliminate the norms of the derivative of the dual solutions zi,Z2- 
In practice, exact solutions 21,^2 are not available. Therefore one way is to use 
accurate numerical approximations for 21, Z2, and further the derivatives can be 
approximated by corresponding difference quotients, see [5] and references therein 
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for practical methods for hyperbolic equations. Another possible way is to eliminate 
these terms by means of stability estimates. 

Example 1. For example, let assume that the desired output functional L*{e) is 
||ei(T)||. One way, see [S] for other alternatives, is to set ji = j2 = zj = and 
zf = ei(T) in (|5.6p . Then, with a = 0, /3 = j = 1, and using the stability estimate 
()2.14p with / = 0, we obtain the a posteriori error estimate 



lei(T)|l<C||lC/i(0)-^.°|l + |l/io(C/2(0)-z;°)|l 

N 



E / {l|f>i-t^2||+||/i„(f>2-./)||+C«(l)( E 



^ 1/2 

hK\\9d-g\\oK 



2 

\dK 



1/2 



+ (c„(i) + c„,^(i))( E '^i'li'-'^i 

+ ||£;fe(f7i -[/2)|| +fc„||-Bfc^(7i|| +fc„||-Bfc / JC{t~s)AhUiis)ds\\ 

Jq 

1/2 , 

dt 



fc„Pfc/ll+fc„( E '^K\\Ek9\\lK) } 



Remark 1. We note that when the convolution kcnel is slightly more regular such 
that ||/C||l2(r+) < oo, following the proof of (|5.15p . we can replace (n N{f3) in (|5.6p 

by 

JC\s^t)ds) E^l^™- 

j=n 

Remark 2. In the a posteriori estimate (|5.6p . with /3 = 2, we have Cn(2) = 1, but 
still Cn.Af(2) can be big when the spatial meshes change. We recall that the error 
representaion (|4.3p has been used for (|5.8p . If we use the error representaion (|4.2p , 
instead, the error indicator 



using the notation ()4.ip . is replaced by 
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Indeed for II5, using the Caucliy-Schwarz inequality and (|5.3[) . we can obtain 

N n „ „t At 



115 = VVV / / IC{t-s)rd{s)ds,Ei,Z2) dt 

N „ n tj/\t 

^Y. E E i / ^(i-s)rd(s) ds \\EnZ2\\ai 
«=i J=i A'eT;;- •'^'-^ 

M n Ti 2 /'2 \ l'2. 

^E / E( E h\m.r,ymiK) ( E ~^k\\e. 

n=l ^"j = l K^f^ Ken 
N „ n 

E / E^™( E ^Aii(^*''^)''wiik) 

n=l"'^", = l f^^-r_J 



ft ^2 II OA- 



< CmaxllV^zaf 

[0,T] 



A-GT 



Hence with (3 = 2 we get the optimal order error indicators, though this error 
indicator is not space-time cellwise. 

Remark 3. We note that for the error estimate ()5.6p there are two types of restric- 
tion on the triangulations; One by Cn.A, that measures the quasiuniformity of the 
family of triangulation, and the other by 5jr^, that measure the regularity of the 
family of triangulations in a slightly different sense. Although maybe not explicitly, 
but Cn,Af (5jr^ can be related. In practice we use finitely many triangulations, 
that means quasiuniformity holds, though possibly with big Qn,Ni see Remark [2] 
This means that we still can use the a posteriori error estimate ()5.6p . But when 
(5jr^ is not sufficiently small, the error estimate (|5.6p does not hold. This calls for 
using local interpolants instead of global L2-pi'ojection Vh- In the next section we 
present an a posteriori error estimate using interpolation, both in space and time, 
on each space-time cell. 



6. A posteriori error estimates based on local projections 

We recall the decomposition of the space-time slab il" = x /„ into cells X" ^ 
K X In, K 7^". Let Ifik be a standard interpolant, for example linear in space and 
constant in time, such that the following error estimates hold for the error operator 
EhkV = [Ihk - 

(6.1) WEhkvU'^ < C(/i^||V'^f IIa'" + kJv\\K^), r = 1, 2, 

(6.2) \\VEhkv\\Kr^ < C{hK\\VMK'^ + kn\\Vv\\K'^)- 

We recall the mesh function hn from (|3.2p . and hmax,n — maxAg-f" hK- We will 
also use the fact that 

(6.3) \\v\\l.= Mt)fdt<KTaa^^\\v{t)f. 

Theorem 3. Let u and U be the solutions of (|2.6p and (|3.4p . respectively, and we 
assume that K- e L2(M+). Then, with e ~ U — u, we have the weighted a posteriori 
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error estimate, for a = 1,2, 

|L*(e)| < Cmax{||V"zi(t)||, ||V2z2(<)||, IliiWIUIisWII, ||Vi2W||} 

(6.4) JX. . 

X |To + 5](T„,i+T„,2)|. 

n=l 

w/iere 

To-||/io"(C/i(0)-^°)|| + ||/ig(C/2(0)-T;")||, 



2)||0"+fc„||f/l-(72|k2. 



+ ll^'(f>2-/)||0" +fc„!|f72-/j|0" 



1 1'^ 

+ kn(^ E ^iflkdllsK") +^"( E ^A'llff/i-ffllaA'") [> 

T„,2- / EA:j/^l|/C|k,(,^.)|( E /^Klk<iWllk')'^' 
Proof. We write the error representation (|4.3p as 

5 TV 5 

(6.5) i*(e)= E 0o,A' + EE E e^^ = io + Ei- 

First we estimate Iq. To this end, recaUing 9o,_r- from (|4.5p . we use the Cauchy- 
Schwarz inequahty and the interpolation error estimate (|6.ip to obtain, for a = 1,2, 

E (C/i(0)-^/°,i?..zi(0))^< E \\Ui{Q) - u''\\K\\Ehkzm\\K 

<c E \\um-A\Kh],\\\7'^z,m\K 

KeT» 

/ --^ \ 1/2 

<q|v"zi(o)||( E h'K\\Ui{Q)^A\K 

Similarly we have 



E {U2iO)-v'',EhkZ2iO))j,<C\\\/^Z2m\\\hl{U2{0) 
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From these two estimates we conclude 

(6.6) lo < Cmax{|lV«zi(0)|l, \\W^Z2m\}{\\h^{Ui{0) - u") \\ + \\hl{U2{0) - «°)||}. 



For the next term, using the Cauchy-Schwarz inequality and the error estimate 
6.ip . we have, for a = 1, 2, 



N 

Il<51 I] \\Ul-U2\\K4EhkZl\\K^ 
N 

E l|f>i-C^2||A-"(/iKl|V"zi||A'"+A:„||ii||A-") 

n=l K&r;: 

" 1/2 / s 1/2 



+ cJ2h E \\ui-u2\\h^) ( E ii^iIIa 

N 

< Cj2{\\KiUi - C/2)|b"|lV"zi|b" + fc„||J7i - [/2|ln"|lii!ln"}, 

n=l 



that using (|6.3p we have 



(6.7) 



Ii <Cmax{|lV«zi(0||J|ii(t)||}Efcy'{||/ij:(f7i-C/2)|lo" 

' n— 1 



In the same way we obtain 



(6.8) 



N 



n—l 
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Now for I3, we use the Cauchy-Schwarz inequality, the scaled trace inequality (|5.5|) . 
and the error estimates (I6.1I)~(I6.2I) to obtain, 



TV 

I3 < X] X] hd- g\\dK^A\EhkZ2\\dK^ 

n=i Ker;: 

N 

<CY, hd- g\\dK'^{hK^^^\\E!,kZ2\\K^^ +hT\\^EhkZ2\\K^^} 
N 

<C!Y.J2 \\9d-g\\aK^ 

X {2/l^/'||V2z2|k" + ^'fcn||i2|k" + hTkn\\yZ2\\K'^} 
^ ( 1/2 

^C-Ei ( E hU9d~g\\lKr^) \\y'z2\\n^ 

n=l ^ KeT^ 



1/2 

+ ^«( E llSd-ffllk'") p2|ln" 

1 /2 

+ fc„(^ E /lArlLgd - glllxn) ||Vi2||f: 



that using (|6.3p we have 
(6.9) 

h<Cmn^{\\V'z2{t)l\\z2{t)l\\Vz2m} 

^ ^ ■ 1/2 / _ s 1/2 



f 1/2 

E^'/'|( E M<\\gd-g\\lKr^) +fc„( E ll.9<i-5iK" 

1 /2 



And similarly 



U<Cma^{\\V'z2{t)l\\z2{t)l\\yz2m} 



r 1/2 



kn[ E ^A'll'-rfllaK' 



1/2 



A'eT" 
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Finally we find an estimate for I5 which includes the convolution terms. First, 
recalling the definition of from ()4.5p . we can write I5 as, 

N T 

15 = -^/ / /C(s-t) ^ (rd{t),EhkZ2{s))gj^dsdt 

Then, using the Cauchy-Schwarz inequality and the assumption /C G L2{R'^), we 
have 

N „ N 



|I5|<^ f Yl W'^'iimaK r IC{s^t)\\EhkZ2{s)\\aKdsdt 

N.N „tj 2/2 

^Y YY \\^d{t)hK( / je{s-t)ds) 



*' \ 1/2 

2 ■ ^ ' 



\\EhkZ2{s)\\%Kds) dt 
^E/ E E \\rm\aKdt\\JCh,^j^^\\E,-,kZ2{s)\\ 

n=l ^'^ i=n i^c-f-i 



N ^ N 

that, using the scaled trace inequality (jS.Sp . the error estimates (j6.ip - (|6.2p . and 
the Cauchy-Schwarz inequality over the triangles, we have 

y '^Z2{s)\\n, 

Ken 
1/2,^ . 



N N , 

|i5i<^E / Eii^iii^a.){( E /^i^ii'-^wiik-) iiv^ 

,1=1 7=ri ^ k-C-rJ 



E '^K'HrdimiK} p2(s)ii 



E 'hKhdimlK) \\^z2{s)\\n.\dt. 
Hence, by (|6.3p . we conclude 



A „ A 

|l5|<Cmax{||V2z2W||,p2(i)IU|Vi2(i)||}E / E II^IU^U.)^. 

' ' ' n=l j=ri 



(6.11) E ^^fii'^rfWiisK) E ^if'ii^'iWi 



Putting estimates (|6.6p - (|6.1ip in (|6.5p we conclude the a posteriori error estimate 
(|6.4p . Now the proof is complete. □ 
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Remark 4. We note that, for example, we can compute llrdljaiC" as 



\rd\\dK^ = 



< 



■rd(tn-l) + 



t — tn-1 



rd{tn)\\Qj^ dt 



1/2 



V3 



[Wrditn^lWaK + WrditnWaK 



1/2 



< 



lkl/^{\\rd{tn-i)\\dK + \\rd{tn)\\dK)- 



Remark 5. Wc recall that to prove the a posteriori error estimate (|6.4|) . we used 
the second error representaion (|4.3p in (j6.5p . Therefore, the error indicator T„ 2 
is not space-time cellwise, since the convolution integral applies to the error term 
EhkZ2- Besides, we need the assumption JC € L2{^'^) for the kernel, that does not 
apply to weakly singular kernels for which we have JC € Li(IR+). Hence for this 
case, when JC E the error indicator Tn,2 is not fulfilled, and we need to 

use the first error representation ()4.2p in (|6.5p . Now we just need to estimate the 
new I5, that using the notation (|4.ip is writen in the form 
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Then, by the Cauchy-Schwarz inequahty and the scaled trace inequahty (|5.5|) . 
we have 



|l5l<CE / E ( E hj,'\\{IC.r,y{tr 

1 /2 

X ( E ^^K\\EhkZ2mlK) \ dt 

^^E/ E ( E K'm-r.ym 



1/2 



1 /2 

E \\EhkZ2{t)\\\ + h],\\\/EhkZ2{t)fK) \dt 

Ken 

^^E/ E ( E K>m-r,ym%K) 

"'max.') 
TV „ n 

CE / P.fc^2WllE( E 



+ / \\\7E^kZ2mJ2~^rna.,,{ E ^^k' II * ''rfF W II k) 

and using the Cauchy-Schwarz inequahty in the integrals over /„, we have 

N / „ / n , ,„x 2 \ 



If/ l/2\ 
I5I < CEll^''fc^2|b" / 

TV 

+ cY,\\yEhkZ2\ 

E^™--^( E /^i^'iK^^'^^^FWi 



n=l 

2 \ 1/2 



lax 



Since by the error estimate (|6.ip we have 



/■ \\EukZ2{t)rdt= I Y. \\EHkZ2{t)\\idt= j2 \\EHkz2mi. 

<C (/*Kl|V'^2||^.+fc^||i2||i.) 



<C{ht,,jV'z2\\l.+kl\\z2\\l.), 
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and similarly by (j6.2p we have 



I \\VEhkZ2{t)fdt= I \\'^EhkZ2{t)\\\dt= \\^EhkZ2m 

■J I. Ji^Ker,: Kef- 

<C (/*^l|V'22|||„+fc^||Vi2 



1^- ) 
if"/ 

<C(/lL..„l|V^Z2||^„+fcfj|Vi2||2„), 

then, recalling (|6.3p . we conclude the estimate 

I5 < Cmax{||V2z2(t)!|, ||i2(i)||, ||Vi2(t)||} Y ^n^'l^^™ + 



?l=l 

2 \ 1/2 

2 

ax 



^ 1.1/2/'/; , , \ 

^ / ^ \"'max.n ~ '^n J 

X ( / (f]Ca.,,( E hKm*r,y{t)\\l,,y^') dt 



1/2, 



This means that, for the case K, E ii(IR+), the error indicator T„,2 in the a poste- 
riori error estimate (|6.4p is 

^ / „ , \ 1/2 

/") .A 

2 \l/2^ 



L \-/^"Vj=i ^g^^. 



2 



1/2 



We note that the last a posteriori error estimate presented in (|6.4p . docs not 
have the restrictions that were mentioned in Remark [3l 



7. Conclusion 

In this work, a space-time continuous Galerkin finite element method has been 
applied to a hyperbolic type integro-differential equation. To provide the main tools 
for adaptive strategies, mainly based on the DWR approach, we have presented 
three error representations in Theorem 1. The main difference between these a 
posteriori error representations is that we apply the convolution integral either on 
the residual or on the error term EhkZ2- Obviously we choose the cheaper one, 
depending on the method we choose for computing or estimating EhkZ2- 

Evaluation of a posteriori error representations require information about the 
dual solution. One way is to find bounds for the dual solution in certain Sobolev 
norms, that can be estimated, e.g., by means of stability estimates. In Theorem 2 
we have presented a weighted global a posteriori error estimate, that is used when 
the goal functional is global. The error indicators have been separated, such that 
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it can be used as a basis for independent adaptation of the time step and spatial 
mesh. However, for this global error estimate, the triangulations are required to 
be quasiuniform and also fulfill a certain kind of regularity, see Remark [3) This 
restriction and also when we use local goal functionals call for local projections. In 
Theorem 3 standard local projections have been used to obtain a posteriori error 
estimates. 

For adaptive strategies the error representations can be used to accurately esti- 
mate the actual errors for getting a stopping criterion of the adaptation process. 
Then we can use the global/local a posteriori error estimates for steering the adap- 
tation process. 

Sparse quadrature can be used to overcome the problem with the growing amount 
of data, in the convolution term, that has to be stored and used in time stepping 
methods. However, we note that this is not an issue for exponentially decaying 
memory kernels, in linear viscoelasticity, that are represented as a Prony series. 
This is due to the existence of a recurrence formula for history updating. We plan 
to address adaptivity strategies based on sparse quadrature together with numerical 
adaptation methods such as the DWR approach, using the theory presented here, 
in future work. 
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